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UNIFORMLY STRESSED INFINITE HOMOGENEOUS MATRIX 
by Alexander Mendelson 
Lewis Research Center 

SUMMARY 

Equations have been derived and a method has been presented for performing an 
elastoplastic analysis of a system consisting of a circular cylindrical inclusion in a homo- 
geneous matrix uniformly stressed at infinity and in a condition of generalized plane strain. 
The material properties of the inclusion and the matrix, including their stress-strain 
curves, are assumed to be arbitrary and independent of each other. Several examples, 
including the limiting cases of a hole and a rigid inclusion, are presented. It is shown 
that the constraints imposed by the rigid inclusion sharply reduce both the stress and 
strain concentration factors over those for the hole. The results for a system roughly 
approximating a graphite fiber in a resin matrix indicate that the fiber acts nearly as a 
rigid inclusion and that strain hardening properties of the matrix play only a minor role 
in determining the plastic strain concentration factor for the case of plane strain. 


INTRODUCTION 

The importance of composite materials in aerospace applications, due primarily to 
their potentially high strength to weight ratios, is well known. The actual design strength 
of such materials is, however, frequently much lower than their potential strength. One 
reason for this is that there may be load components perpendicular to the fiber axes. The 
fibers then act as inclusions producing stress concentrations in the matrix, which is 
already considerably weaker than the fibers. 

An elastic analysis on the effect of a cylindrical inclusion in a homogeneous matrix 
is presented in reference 1. Herein a more realistic elastoplastic analysis is made to 
determine the stress and strain concentrations due to such inclusions. A complete anal- 
ysis would have to include a large number of irregularly shaped and spaced inclusions as 
well as the orthotropic character of the fibers. The present analysis is therefore merely 



a first step, wherein a single circular inclusion in an infinite homogeneous matrix, uni- 
formly loaded at infinity, is considered. 

The method used in the analysis is the successive approximation or iterative tech- 
nique for solving elastoplastic problems as outlined in references 2 and 3. This tech- 
nique has been successfully used for the problem of an infinite plate with a hole in refer- 
ences 4 and 6 and for a rigid inclusion in reference 7. In both cases, the plane stress 
problem of a thin plate was considered. In the present analysis, the generalized plane 
strain problem of an infinitely thick plate is treated for arbitrary matrix and inclusion 
properties. Furthermore, Hencky's total plasticity theory is used. Although the incre- 
mental theory can be used with equal ease, the results shown in reference 6 indicate that 
for this type of problem there is negligible difference between the total and incremental 
theory. By using total plasticity theory an appreciable amount of computor time can be 
saved. 

The methods used in this report were developed by the author while at NASA Lewis 
Research Center. The bulk of the research was carried out while the author was on leave 
from Lewis at Case-Western Reserve University. The author acknowleges the sponsor- 
ship during this period of the Advanced Research Projects Agency, Department of Defense, 
through a contract administered by the Air Force Materials Laboratory. This report is 
also being issued as an Air Force Materials Laboratory Technical Report. 


SYMBOLS 


a 

A, B 

E 

e 

KpK 2 , . . - Kg 

m 

P,Q,R 

r 

S 

S.O 

T 

U 


radius of inclusion 
integration constants 
elastic modulus 
ratio of strain to yield strain 
constants defined by eq. (18) 

ratio of slope of linear strain-hardening curve to elastic modulus 

plastic strain functions defined by eq. (4) 

radial coordinate 

ratio of stress to yield stress 

ratio of applied uniform stress at infinity to yield stress 
ratio of average axial stress due to end loads to yield stress 
dimensionless radial displacement defined by eq. (1) 
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u radial displacement 

a ratio of matrix yield stress to inclusion yield stress 

j3 ratio of matrix elastic modulus to inclusion elastic modulus 

e strain 

p, Poisson's ratio 

p dimensionless radius, r/a 

p max maximum value of p used in computations 

p radius to plastic zone boundary 

r 

cr stress 

Subscripts: 
e equivalent 

I inclusion 

M matrix 

o yield 

r radial 

z axial 

9 tangential 

Superscript: 
p plastic 


ANALYSIS 

General Relations 


Consider an infinite plate uniformly loaded at infinity as shown in figure 1. Be- 
cause of the symmetry, the problem is most simply cast into polar coordinates. 




Figure 1. - Infinite body with circular inclusion. 


The dimensionless quantities used in this report are defined as follows: 



where o Q and e Q are the yield stress and yield strain, respectively (these will be more 
precisely defined later), a is the radius of the inclusion, and u is the radial displace- 
ment. The equilibrium, compatibility, and stress-strain relations can be written as 
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( 2 ) 


dS S a - S ^ 

r _ 6 r 

dp p 

de g _ e r - e 6 

dp P y 

e r = S r “ M(S 0 + S z ) + eP 

e e = s e- ^ (s r + s z ) + e 0 

e = S - p(S + S fl ) + = constant 


Equations (2) can be integrated to give the following relations: 


S_ = A + — + 


2 9/i 2, 

p 2(1 - p ) 


1 P( p ) - J— -?£_ Q(p) 


2(1 - n ) 


S q = A - — + I 

0 2 


P(p) + -I ' Q(p) 


R(p) 


P* 2(1 - p‘) 2(1 - p*) 1 - p‘ 

S z = p( S r + S 0> + e ? + e 0 + e z 


e r = (1 - M - 2p 


)A + (1 + p) A + I 1 -^fp(p) + Q(p)l + eP - pe 

p 2 2 1 - p L J i - p r 

0 = (1 - P - 2p 2 )A - (1 + p) + ILlIeVp ( P ) + Q(pjl - pe 

2 2 1 - p L - 1 

P 

/*p /*p 

+ s e^ p dp + J Q ( e ? + e e)P d P 

u = (1 - p - 2p 2 )Ap - (1 + p) ® + 1 p [p (p) + Q(p)] - ppe z 


e„ = T - 2 lim — 
z 2 
P-°° P L 


( 3 ) 
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where 


P(p) s 







R(p) = M e ? - (! " M) e 0 


( 4 ) 


and T is the dimensionless average axial stress due to the end loads. The constants 
A and B will depend on the boundary conditions as will be shown, and the lower limit on 
the integrals c will equal zero for the inclusion and one for the matrix. 

The plastic strains are related to the total strains by the modified Prandtl-Reuss 
relations (ref. 2). 


A 


e r = ~- ( 2e r " e n ~ e J 


3e„ 


V=~ ( 2e fl - e r - O 


3e 


(5) 


where 


e 


e 



,2 

q) + 



+ ( e Q 



( 6 ) 


and e is related to e through the stress-strain curve and the relation 
D ^ 


e D = e - - (i + n)S 
p e 3 e 


(7) 


where S , the dimensionless equivalent stress, is a function of the equivalent plastic 
strain e p and is the ordinate of the dimensionless uniaxial stress-strain curve as shown 
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V e p + 3 (1 + |1,S e 


(a) For nonlinear strain hardening. 

/ 

/ 

/ 

/ 



1 


Ratio of strain to yield strain, e 


(b) For linear strain hardening. 

Figure 2. - Dimensionless uniaxial stress-strain curve. 


in figure 2(a). 
to give 


For the case of linear strain hardening, equation (7) can be solved for e 


6 P = 


e e - - (1 + p) 
e 3 

1 + (1 + p) 

1 - m 3 


( 8 ) 


where m is the ratio of the slope of the stress-strain curve in the plastic range to the 
elastic modulus, as shown in figure 2(b). For a perfectly plastic material, equations (7) 
and (8) reduce to 


e = e - - (1 + p) (8a) 

p e 3 

The above equations are valid both in the inclusion and in the matrix. The constants 
A and B will be different for the matrix and for the inclusion as will the material prop- 
erties. We will henceforth distinguish between the matrix and the inclusion by using the 
subscripts M and I for matrix and inclusion, respectively. 


Evaluation of Constants 

An evaluation of the constants makes use of the boundary conditions. As previously 
noted for the inclusion 


c = 0, 0 < p < 1 


and for the matrix 


c = 1, 1 < p 

It also follows from the first or second equation (eq. (3)) that, for finite stresses at 
P = 0, 


B 


I 


= 0 


(9) 


Substituting the values of S r and S 0 from equations (3) into the expression for e z 
and taking the indicated limit gives 
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e = T - 2n. A,. - 
z M 


M 


M 


i 2 

1 ■ m m 


P M<”> 


Let the dimensionless stress at p = °° be designated by 



°oM 


(10) 


where a QM is the yield stress of the matrix. Then, from the first of equations (3), 

1 2 

A M = S ~- ~ 2 P M (0 °) + ~~ 

^ - 4> 2d - 4> 

where, from equations (4), 


Hence, 


/ °° P P 

erM ~ e0M dp 

q m (°°) = 0 


( 11 ) 


A M = S -c P M<"> < 12 > 

2 d - 4> 


Note that, since the plastic strains die out as p increases, the upper limit in equa- 
tion (11) may be replaced by some arbitrary radius p where p is greater than 

max max 

the radius of the plastic zone. Thus, 


p mH 



p 



dp 


(13) 


p m ax — radius of plastic zone 

The two remaining constants Aj and B M can be determined from the conditions 
that the radial stress and displacement should be continuous across the matrix- inclusion 
interface, that is, 
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CT rI< a > = CT rM (a ) 
V a) - u M (a) 


s rI (i)=!^s rM (i) 

CT oI 


s r i(D = «s rM (i) d' 

where 

CT oI 

is the ratio of the matrix yield stress to the inclusion yield stress. Note that S rI is 
defined as the ratio of the inclusion radial stress to the inclusion yield stress and S rM 
is defined as the ratio of the matrix radial stress to the matrix yield stress. Similarly, 


/SUjd) = «u M (i) 


(15) 


where 


„, E m „ _“l E I „ _ u M E m 

E x aa Ql aa. 


(16) 


'oM 


and E and Ej are the elastic moduli of the matrix and inclusion, respectively. 
From equations (3), (9), (12), (14), and (15) and noting that Pj^(l) - - 

Aj = Kfoo + ^(D + + K 4 P M ( °°> + K 8 e z ^1 

b M = K 5 S oo + + K 7 P M (o0 > + ^ K 8 e z 


(17) 


where 
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K, 


1 + M M + (! + M x )( 1 " 2 Pj)i3 

2(1 - nf) 

1 - 2/ij /l + m m 




- 0 


K 3 = 


K 4 


2(1 “ Mj) y 1 + Mj 

1 + Mm + (! + Mj)(l - 2 Mj)0 


-a 


1 + M M + (! + MpC 1 “ 2 V I )P 
K 1 

k 5 =— - 1 

a 




K 6 = 


(1 - 2 jLLj)/3 


K, 


a 


K c 


K 7 = - 


2d - Mi) 


K 8 = (^M-^l) K 4 


(18) 


Note that the coefficients through Kg are functions of the four material constants 

Mjj Mjyp a an< ^ 0* 


Computation Procedure 

The stresses, total strains, and plastic strains can now be computed by the succes- 
sive approximation or iterative method described in reference 2. The inclusion and the 
matrix are each divided into finite radial intervals. For the inclusion, the stations range 
from zero to 1. For the matrix, the stations range from p = 1 to p = p max where 
Pmax I s an arbitrary dimensionless radius larger than the radius of the plastic zone. 
Initially, the radius of the plastic zone is not known; therefore, «one must guess at a 
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reasonable value for p . It may be desirable to change this value after the first iter- 
ation so that most of the stations are in the plastic zone. 

The coefficients through Kg are computed once from equation (18). Assuming 
all the plastic strains to be zero so that Pj, Qj, P M , and Q M are all zero, Aj and 
B m are computed from equation (17), A M from equation (12), and e z from equation 
(10). Equations (3) now give the complete elastic solution for the inclusion and the ma- 
trix. Beginning with the elastic solution, the iterative scheme for obtaining the elasto- 
plastic solution proceeds as follows: 

(1) Using the values of total strains just computed, the equivalent total strain e g is 
calculated at every station of the inclusion and the matrix by equation (6). 

(2) The equivalent plastic strain ep at every station is determined from equation (7) 
and the stress-strain curve (or eq. (8) for linear strain hardening). If ep at any station 
is less than zero, there is no plastic flow at that station and the plastic strains are set 
equal to zero at that station. 

(3) The individual plastic strains at every station are computed by means of equa- 
tions (5). 

(4) The plastic strain integrals Pj, Qj, P M , and Q M are computed from equa- 
tion (4) for every station. 

(5) e z , Ajyp Aj, and B M are calculated from equations (10), (12), and (17). 

(6) The stresses and strains are then computed from equations (3). 

(7) Return to step (1) and continue iterating until convergence is obtained, that is, 
until two successive solutions differ by less than some arbitrarily preassigned value. 

The outlined computation scheme, which was programmed for a digital computer, 
gives rapidly and accurately the complete stress and strain fields both in the inclusion 
and the matrix. Both the elastic properties and stress-strain curves will, of course, 
be different in the two media. 


RESULTS AND DISCUSSION 

The technique presented will now be illustrated for three different cases: The two 
extreme cases of a plate with a void and a plate with a rigid inclusion, and the case of a 
relatively weak matrix with a relatively strong inclusion. In all cases, the condition of 
plane strain with e z = 0 is assumed. 


Plate With Hole or Void 


A plate with a hole represents the limiting case 


a — °o 


oo 


( 19 ) 


P - 


The calculation proceeds as described except that Aj need not be computed. The stress 
and strain fields are obtained only in the matrix. Some results are shown in figure 3. 

In figures 3(a) and (b) the effective stress and strain concentration factors are 
plotted against the applied stress at infinity for various values of the strain hardening 
parameter m^.. The effective stress concentration factor is defined as the ratio of 
effective or equivalent stress at the edge of the hole to the equivalent stress at infinity; 



Figure 3. - Effect of strain-hardening parameter on plane-strain case of infinite plate with 
hole. Ratio of matrix yield stress to inclusion yield stress, «>; ratio of matrix elastic 
modulus to inclusion elastic modulus, 
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Dimensionless load, 


(b) Effective strain concentration factor. 
Figured - Continued. 


that is, 

K - S e (1) _ V 1 * 

CT s e(°°) - V- 


( 20 ) 


where 


S e = 


(s r - SJ* + (Sg - S ; + (S - S )‘ 


1/2 


The elastic stress concentration factor as defined is readily shown as 


^(elastic) = . / l + 

y a - 2 h m > 2 

For a value of ju M of 0. 35, which is the value used in these calculations, the elastic 
stress concentration factor is equal to 5. 86, which is the common starting point of all 
the curves in figure 3(a). 

Note that, for the case of plane stress, the stress concentration factor is independent 
of Poisson's ratio and is equal to 2 for the elastic case. The large initial value of the 
stress concentration factor for the case of plane strain does not mean, however, that 
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plastic flow will start at a lower load than for plane stress. On the contrary, it can 
readily be shown that the equivalent stress S g will equal 1 at the hole for a value of 
= 0. 569 compared with = 0. 5 for plane stress. For a plate without a hole, 
must equal 3. 33 (for p = 0. 35) in order that S g = 1, whereas, for the plane stress case, 
S g = 1 when 8^=1. Thus, it requires more than 3 times the load for yielding to occur 
in the plane-strain case compared with the plane stress case. However, the effect of a 
hole is felt much more for plane strain than plane stress, so that the hole weakens the 
plate by a factor of almost 6 for plane strain compared with a factor of 2 for plane stress. 
The result is that the load to produce yielding at the hole is equal to 0. 569 for plane 
strain compared to 0. 5 for plane stress as indicated previously. The plane-strain case 
still requires a larger load for yielding to occur than the plane-stress case, but the dif- 
ference is no longer as great when there is a hole present. 

Whether definition (20) is appropriate depends on how the stress and strain concen- 
tration factors are to be used. The concentration factors, by definition, give the ratios 
of the stress or strain with the hole present to the corresponding values without the hole. 
They therefore represent the weakening effect of the hole. The equivalent stress and 
equivalent strain are used herein to represent the concentration factors because they 
are the determining quantities in plastic flow. If some other quantity is more useful in 
determing failure, then that quantity should be used rather than the equivalent stress or 
equivalent strain in defining the concentration factors. 

Figures 3(a) and (b) show how the stress concentration factor decreases with load 
and the strain concentration factor increases with load because the material is flowing 
plastically. The importance of strain hardening is readily evident. It is interesting to 
note that, although without strain hardening the stress concentration factor decreases 
monotonically with load, the presence of strain hardening introduces a minimum point 
in the stress concentration as function of load curve. For a high enough load, the strain 
hardening of the material is sufficient to overcome the stress relaxation due to plastic 
flow, and the stress concentration factor starts rising. This effect does not occur for 
the plane stress case (see ref. 6), probably, because the additional axial constraint 
inherent in the plane-strain case is missing. 

The strain concentration factor used in figure 3(b) is defined by 


K 


e 


e e < " ) |(1 + I * M )<1 - 2 Mm )S 


( 21 ) 
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(c) Plastic zone radius. 


Figure 3. - Concluded. 


For the elastic case, 


e 0 (l) =- (1 + m m ) 


3 + (1 - 2 Mm ) 


1/2 


so that, for = 0. 35, 


K £ = 5. 86 = K a 


Figure 3(c) shows the growth of the plastic zone with load. Again the effect of strain 
hardening is very evident. The plate becomes completely plastic when S M equals 3-|. 


Plate With Rigid Inclusion 


For a rigid inclusion 



a = 1 3 = 0 


( 22 ) 


The results using these values are shown in figure 4. Note the very large difference be- 
tween these results and those for the plate with the hole. The restraints produced by the 
inclusion greatly lower both the stress and strain concentration factors. Furthermore, 
the effect of strain hardening is no longer as great. The plastic zone size is greatly re- 
duced as compared with the hole (figs. 3(c) and 4(c)). 




(b) Effective strain concentration factor. 


Figure 4. - Effect of strain-hardening parameter on inclusion. Ratio of matrix yield stress 
to inclusion yield stress, 0; ratio of matrix elastic modulus to inclusion elastic modulus, 
0 . 
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Strain concentration Stress concentration factor, 

factor, K € 




Dimensionless load, S^ 


(b) Effective strain concentration factor. 

Figure 5. - Effect of strain-hardening parameter on plane-strain case of infinite plate with 
fiber inclusion. Ratio of matrix yield stress to inclusion yield stress, 0. 1; ratio of matrix 
elastic modulus to inclusion elastic modulus, 0.02. 


Plastic zone radius, pp 



Figure 5. - Concluded. 
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Inclusion With a = 0.1, (3=0.02 


As a final example, a matrix inclusion combination with a = 0. 1 and /3 = 0. 02 was 
considered. These properties correspond roughly to a graphite fiber in a resin matrix. 
The results are shown in figure 5. Note the similarity of these results to those of fig- 
ure 4 for the rigid inclusion. Of particular interest is the fact that the strain-hardening 
properties of the matrix are relatively unimportant. Only one curve was drawn in fig- 
ure 5(c) because the difference between the different strain-hardening parameters is 
very small. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 28, 1967, 

129-03-08-04-22. 
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